Optomechanical quantum information processing with photons and phonons 
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We describe how strong resonant interactions in multimode optomechanical systems can be used 
to induce controlled nonlinear couplings between single photons and phonons. Combined with linear 
mapping schemes between photons and phonons, these techniques provide a universal building block 
for various classical and quantum information processing applications. Our approach is especially 
suited for nano-optomechanical devices, where strong optomechanical coupling on a single photon 
level is within experimental reach. 
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Optomechanics describes the radiation pressure inter- 
action between an optical cavity mode and the motion 
of a macroscopic mechanical object, as it appears, for 
example, in a Fabry-Perot cavity with a moveable mir- 
ror [T]. First demonstrations of optomechanical (OM) 
laser cooling [2] have recently attracted significant inter- 
est and led to tremendous progress in the development of 
new fabrication methods and experimental techniques for 
controlling OM interactions at the quantum level. Apart 
from ground state cooling [3, i4] , this includes the demon- 
stration of slow light 16] , and the coherent interconver- 
sion of optical and mechanical excitations [7l Hj . These 
achievements pave the way for a new type of quantum 
light-matter interface and give rise to interesting per- 
spectives for novel OM based quantum technologies. As a 
solid state approach, such an all OM platform would ben- 
efit directly from advanced nanofabrication and scalable 
integrated photonic circuit techniques. In particular, as 
compared to atomic cavity QED [5] setups, no trapping 
of particles (close to surfaces) is required. 

In this work we study strong OM coupling effects in 
multimode OM systems (OMS) and describe how reso- 
nant or near-resonant interactions in this setting allow 
us to exploit the intrinsic nonlinearity of radiation pres- 
sure in an optimal way. Our approach is based on the 
resonant exchange of photons between two optical modes 
mediated by a single phonon. This resonance induces 
much stronger nonlinearities than achievable in single 
mode OMS, where due to the frequency mismatch be- 
tween optical and mechanical degrees of freedom only 
off- resonant couplings occur [TUl - fT5| . Consequently, mul- 
timode OMS provide a promising route for accessing the 
single photon strong coupling regime [12] , where the cou- 
pling f^o well as the mechanical frequency exceeds 
the cavity decay rate k. Using state of the art nanoscale 
devices [H [T3] , this regime is within experimental reach, 
and here we discuss how strong OM interactions in a 
multimode setup can be used to generate single pho- 
tons and to perform controlled gate operations between 
photonic or mechanical qubits. Combined with very re- 




FIG. 1. (color online) a) Setup of two tunnel-coupled OM 
crystal cavities (see Ref. (4] [14] for more details), b) Level di- 
agram showing the lowest mechanical and optical excitations 
in a two mode OMS. Resonant coupling {5 = 0) occurs when 
the tunnel splitting 2J between the optical modes is compa- 
rable to the mechanical frequency ujm- c) Different sets of 
strongly and weakly coupled optical modes and control laser 
fields can be used for nonlinear interactions (tJs , i^a , i^l ) and 
purely linear photon storage and retrieval operations {uj'c,u)'j^). 



cently developed photon-phonon interfaces and quantum 
memories based on linearized OM couplings [71 [H [15], 
our results provide a basis for efficient OM classical and 
quantum information processing with applications rang- 
ing from photon transistors to quantum repeaters and 
networks. 

Model. We consider a setup of two tunnel-coupled 
OMS [TSHTO] as schematically shown in Fig. [l] focusing 
on the OM crystal design [H [T3| as a specific example. 
Each OMS i = 1,2 is represented by an optical mode of 
frequency Wc and a bosonic operator Ci , which is coupled 
via optical gradient forces to the motion of an isolated 
mechanical resonator bi with vibrational frequency w^. 
The total Hamiltonian for this system is [fi — 1) 

H ^in4bt + ^cc\ci + goc\ci[bi + 6j) 

(1) 

- J{c\c2 + cic\) + ^ r!Kc,e'"^* + H.c), 
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FIG. 2. (color online) a) Energy level diagram of a resonantly 
coupled OMS, 5 = 2J — ujm = 0, and for a single mechan- 
ical mode in the ground state, b) Excitation spectrum and 
(?'^''(0) for a weak coherent field exciting the Ca mode, where 
<?o/k ~ 8 and no — Q^^jK^ . c) Minimal value of (;'^-'(0) as a 
function of the OM coupling strength go and different values 
of Nth- The analytical results (solid lines) given in the text 
are in good agreement with exact numerics (circles). The 
dashed line shows the asymptotic scaling ~ 8k^/(?o at zero 
temperature. 



where J is the tunneling amplitude between the opti- 
cal modes and go denotes the single photon OM cou- 
pling. The are the local amplitudes of an exter- 
nal control laser field of frequency a;^. Below we also 
consider an additional set of cavity modes and driving 
fields with frequencies and w^, respectively. As in- 
dicated in Fig. [ijc) we assume these modes to be sepa- 
rated in frequency and used for cooling the mechanical 
mode [201 m], and linear photon storage and retrieval 
operations [3 H [H [23] only. 

Apart from the coherent dynamics described by 
Eq. ([T]), we include dissipation through cavity decay and 
mechanical damping and model the evolution of the sys- 
tem density operator p by a master equation (ME) 



(2) 



where V\c\p = 2cpc^ - {c'fc,p}+, and £^ = Yl,i ^{Nth + 
l)'D[bi] + ^Nth'Dlbl]. Here k is the optical field decay 
rate, 7 = utm/Q the mechanical damping rate for a qual- 
ity factor Q and Nth = (e'^""*/'^^"^ ~ the mechanical 
equilibrium occupation number for temperature T. Be- 
low we identify = ^{3Nth + ^) as the characteristic 
decoherence rate for mechanical qubit states [23] . 

Resonant strong coupling optomechanics. In the fol- 
lowing we are interested in the strong coupling regime 
'^miSo ^ KjT^, and our main goal is to show how the 
multi-mode OMS described by Eq. ([!]) can be used for 
implementing controlled interactions between qubits en- 
coded in photonic or phononic degrees of freedom. To 
illustrate this we first consider a single mechanical res- 



onator, b = bi, uj„ 



We introduce symmetric and 



antisymmetric optical modes Cs^a = (ci ± C2) /v^ with 
eigenfrequencies Us.a split by 2J. Further, we assume 



that ojm ^ 2J ^ go,K,\5\, where 5 = 2J — uj^ (see 
Fig.jljb)). This condition can be achieved in nano-scale 
OMS where a;„ - GHz [Hini EH [26] and allows us to 
make a rotating wave approximation with respect to the 
large frequency scale ujm 2J. After changing into a 
frame rotating with ujl we obtain |17| 



(3) 



where A^^a = i^l — '-^s,a are the detunings of the driving 
field from the and mode respectively, and H^(t) — 
'^r)=s a i^vi^)^v + H.c.) accouuts for the external driving 
fields with slowly varying amplitudes Qs,ait) = (fliit) ± 

n2{t))/V2. 

The two-mode OM coupling in Eq. ([s]) describes pho- 
ton transitions between the energetically higher mode Cq 
to the lower mode c,, while simultaneously absorbing or 
emitting a phonon. For (A^ — Aa — i^m) = ^ = 0, this 
leads to a resonant interaction between states jria, n,, nm) 
and \na ~ l,ns + 1, Um + 1), where n^, Us and Um label 
the occupation numbers of the two optical modes and the 
mechanical mode, respectively. In analogy to atomic cav- 
ity QED [9], the nonlinear scaling of the corresponding 
transition amplitudes ^ y^na{ns + !)(«„ -I- 1) results in 
an anharmonic level diagram as shown in Fig.[2|^a). If go 
exceeds the cavity linewidth k, one and two photon tran- 
sitions can be spectrally resolved, indicating the onset of 
strong single photon nonlinearities. 

An OM single photon source. As a first application of 
the nonlinear OM interaction we discuss the use of the 
OMS as a single photon source, which is characterized 
by a vanishing equal time two-photon correlation func- 
tion g^^)(0). In Fig. we plot the excitation spec- 
trum {clca) and g'-'^\0) = {clclcaCa) / {clca)"^ , for the 
case where only the Ca mode is weakly driven. Around 
the single photon resonances Aq — ±go/2 we observe 
strong anti-bunching g^^^(O) < 1 as a clear signature of 
non-classical photon statistics. To quantify this effect we 
assume that Tm ^ which allows us to treat subspaces 
connected to different \nm) separately. For weak driving 
fields <C K the system dynamics can then be restricted 
to the six states |0a, Os, n™), 0,, |0q, 1,, -I- 
1), \2a,0s,nm), \laAs,nm + ^), I Oa , 2, , n„ 2) and we cal- 
culate the relevant occupation probabilities pi,o,ra„ and 
P2,o,n„ to leading order in fla [IZ]- We obtain 



(X„(2X„-52)) 

(4) 

where d — Aa ~ in and ^ d^ — g'^{n + 1) . By taking 
the appropriate thermal averages, (tIq) — CnPi.o,™ 

and g^^HO) = 2 Cr^P2,o,n/(?^Q)^, where Cn = (1 - 
^-phu,^y-f}hu,^n p-1 ^ ^^y^ ^j^g ^-^Q photon corre- 
lation function can be evaluated for arbitrary tempera- 
tures T. 



FIG. 3. (color online) A single phonon - single photon tran- 
sistor, a) An incoming photon in port (A) passes through the 
interferometric setup and leaves through port (A) or (B), de- 
pending on the phase shift A<j!> acquired upon reflection from 
the two mode OMS. b), c) For a mechanical system in state 
|0m), the OMS exhibits a single resonance at ujs (A^ — tt), 
while for state \lm) the resonance splits hy go ^ n and the 
photon does not enter the cavity (A^ — 0). 



In Fig. ^c) we plot the minimal value of g'--'^^0) as 
a function of the coupling strength Qq and for different 
Nth- As the OM coupling increases we find that for 
T — the minimum of the correlation functions scales 
as miUA, Ag'-'^H^)} - 8^2/5^. This demonstrates an im- 
proved scaling over off-resonant photon blockade effects 
occurring in single mode OMS [12^. Since the positions 
of the single and two-photon resonances depend explic- 
itly on the mechanical state \nm), finite temperature de- 
grades the quality of the single photon source. Neverthe- 
less, with increasing coupling strength the anti-bunching 
effect is surprisingly robust and can be observed even for 
a few tens of thermal mechanical quanta. 

Single-phonon single-photon transistor. Given the 
ability to generate single photons, Fig. [3] illustrates a 
basic scheme for using the same resonant OMS to im- 
plement a two qubit gate [28^. First, we assume that 
the state of a control photon is mapped onto a mechan- 
ical superposition state ajOm) + /3|ljn). This can be 
achieved with conventional cooling followed by photon- 
phonon conversion techniques using linearized OM inter- 
actions with an auxiliary mode lo'^ (see Fig. [ijc). Next, 
a single target photon of central frequency ^ lOs is sent 
through the interferometric setup as described in Fig. [3] 
If the mechanical mode is in the state |0m) the incom- 
ing photon couples to a single resonant state |0a, ls,Om) 
(see Fig. [T|b) , such that it enters the cavity and picks up 
a phase before being reflected. Instead, if the mechani- 
cal resonator is in the state the resonant coupling 
between |0a. Is, Im) and |la, 0^, Om) splits the cavity res- 
onance and for go > k the photon is reflected without a 
phase shift. Under ideal conditions the final result is an 
entangled state 



a|0™,lA,Os) +/3|1™,0a,1b), 



(5) 



FIG. 4. (color online) a) OM quantum memory, where 'path- 
encoded' photonic qubits are stored in long-lived mechanical 
states using tunable linearized OM interactions ~ i^i{t). De- 
terministic gate operation between stationary qubits are im- 
plemented by a controlled phonon-phonon interaction ~ fig (t) 
as described in the text, b) The total error eg for imple- 
menting a controlled phase gate between two phononic qubits 
is minimized with respect to As and plotted as a function 
of K and Tm (see text). The parameters for this plot are 
go/{2n) = 50 MHz, 7/(27r) = 4 kHz, a = 1 and go/ 5 = 1/3. 



By assuming that the storage and the retrieval of the 
control photon can be achieved with high fidelity, the 
error for producing the entangled state ([5| with a — (3 — 
I/V2 is approximately given by 



(TpK)2 



(6) 



where A and B are the two ports of the interferometer. 
This state can be converted back into an entangled state 
between the initial control and the target photon. 



where Tp is the duration of the single photon pulse. The 
individual contributions in Eq. ^ arise from an imper- 
fect photon reflection, the finite spectral width of the 
photon pulse, and mechanical decoherence, respectively. 
The minimal error is achieved for « \/T?T\^ where 

we obtain e w niax{4:K^ /g^, ■\/^m/^'^}- Assuming an OM 
crystal device with a;m/(27r) = 4 GHz and Q ~ 10^ as 
dicussed in Ref. U, but with an improved OM coupling 
(/o/(27r) = 50 MHz and a lower decay rate k/{2tt) = 5 
MHz, we obtain gate errors e w 0.1 for environmental 
temperatures around T « 100 mK. 

Phonon-phonon interactions. Finally, we consider the 
possibility to perform a controlled gate operation be- 
tween two qubits stored in long-lived mechanical modes. 
Our approach is depicted in Fig. ^a.) , and combines the 
long coherence times of OM quantum memory [71 [51 [221 
I23j with the practical utility of exploiting interactions be- 
tween stationary phononic qubits. We focus on the limit 
r,„ <^ K, and assume that optical (e.g. 'path encoded') 
qubits are first mapped onto long-lived states \0m) and 
|l„i) of two or more mechanical modes. The OM cou- 
pling is then employed to generate nonlinear interactions 
between the phonons only. 

We consider nonlinear interactions between two me- 
chanical modes bi and 62 described by Eq. (!]), detuned 
from resonance such that go < |(2J — w^) and direct 
transitions between photons and phonons are suppressed. 
To obtain the effective phonon-phonon interactions, we 
first diagonalize H to second order in = go/{2J — w^) 



4 



with the transformation H e^^He , where S = 
f (4ca(a&I - 64) - H.c). This yields H = Ho + Hg + 
Hn{t), where Hq = -A^cjcs - Aa4ca + J2t ^U^lh, and 

^9 - f [(4cs+i)4c„(a+6) + (4ca-4cs)A/;] . (7) 

We neglected small corrections to the driving Hamilto- 
nian The phonon operator in Eq. ([t]) is given by 

Mt, = iib\b, + i2b\b2 - (^1 + 6)(4^2 + 46i)/2. For sim- 
plicity we focus on symmetric detuning, uj}-^ ~ 2 J =F J, 
where Mb = ^{b\bi — 62^2)- The transformation also 
modifies the dissipative terms in the Eq. ([2| ; most impor- 
tantly we find an optically-induced decay channel for the 
mechanical modes, + ng^ / {A5'^)'D[cs{bi + 62)]. 

We assume that only the Cg mode is weakly driven 
by a slowly- varying control field ^s{t)- In this case the 
Ca mode remains unpopulated and we neglect it. Next 
we shift the driven mode, — > a -I- Cj by the classical 
amplitude a, yielding an effective ME for c^, 61 and 62- 
Finally, we adiabatically eliminate the Cs mode, valid in 
the limit \a\ ~ 0(l) and {gl\a\lA5) < \l\s+iK\, to obtain 
an effective master equation for the mechanical modes, 

Pm = - i[Hm + J^{b\bl - 4^2)^Pm] + -C^Pm 

i 

Here 7' = K|ap(7§/(2(5^) and the phonon-phonon interac- 
tion and the phonon dephasing rate are given by 

16<52(A2-f ^2)' * 1652(A2 + k2)- 

The effective Hamiltonian in Eq. ^ describes a phonon 
nonlinearity with a tunable strength K{t) ^ \a{t)\'^ . The 
nonlinear interaction is given by 

i/i„t ^ 2Kb%b\b2, (10) 

and when acting for a time tg = 7r/(2A) this Hamil- 
tonian implements a controlled-phase gate between two 
qubits encoded in states \Qm) and During this time, 

phonons experience intrinsic and optically-induced deco- 
herence as seen in Eq. ^ . In Fig. |4] we plot the re- 
sulting gate error = 1 — {4'o\Pmitg)\iJJo) for an initial 
state I^Ao) = idOm) + llm))*^^ optimized with respect 
to As. Using the total decoherence rate of this state, 
Fdccoh = SFm -I- F^ + 772, we find that eg oc Fdccoh/A 
is minimized for jA^I ~ 30/2 where eg oc 4(«;/(7o)- While 
this scaling with 170 is weaker than for a gate based on 
photon reflection (see Eq. ([6|), the ability to perform a 
gate between stationary qubits represents an important 
advantage of this approach. 

Conclusions. We have described single photon and sin- 
gle phonon nonlinear effects in strongly coupled multi- 
mode OMS. We have shown how induced nonlinearities 



on or near resonance can be used for controlled quan- 
tum gate operations between flying optical or stationary 
phononic qubits. Our results provide a realistic route 
towards the quantum nonlinear regime of OMS, and a 
framework for future OM information processing appli- 
cations. 
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Note added. During the final stage of this project a 
related work by M. Ludwig et al. appeared |29) . 
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